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ABSTRACT 

Aims. Astrophysical objects frequently exhibit some irregularities or complex behaviour in their light curves. We focus primarily on hot stars, 
where both radial and non-radial pulsations are observed. One of the primary research goals is to determine physical parameters of stellar 
pulsations by analyzing their light curves or spectra, focusing on periodic or quasiperiodic behaviour. 

Methods. We analyse the feasibility of classical methods for period searches in a nonlinear chaotic system, such as the Rossler system, where 
a period does not exist at all. As an astrophysical application of the chaotic system, we utilize a simple model of stellar pulsation with two 
different sets of parameters corresponding to periodic and chaotic behaviour. For both models we create a synthetic signal, and then apply 
widely used methods for period finding, such as the phase dispersion method and periodograms. For comparison, a quasi-periodic signal is 
employed as well. 

Results. The period analysis indicates periods even for the chaotic signal. Such periods are apparently spurious. This implies that it is very 
problematic to distinguish chaotic and quasiperiodic process by such an analysis only. 
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1. Introduction 

Variability of astronomical sources like stars is an important 
aspect of their observation. We can detect variability in their 
light and/or spectra on different timescales. Dynamical pulsa- 
tions, radial or non-radial, may be one of the main sources of 
intrinsic variabilit y in many types of stars. As was shown by 
lEddingtonl ( ll919l) . stars may exhibit pulsation due to the unsta- 
ble layers in the stellar envelope. The first theoretical modelling 
' of pulsations predicted regular oscillations and such modelling 
' was succesfully applied to objects like classical Cepheids. But 
' later, observations indicated other types of pulsating stars, e.g. 
Population II Cepheids, which in some cases exhibit irregular 
patterns in observational data. 

There was no reliable explanation of such variability un- 



precise one-zone model also produced chaotic solutions, with a 
strange attractor similar to the Rossler one (see lRossler ■ Il976l) . 



Detailed, more realistic calculations of the state equation for 
he lium ionization in a one zone model, which was performed 
by lBuchler & RegevI (Il982h . confirmed that pulsation instabil- 
ity may lead to chaotic behaviour 

Goupil et al. 1 (Il986h studied light curves of the pulsating 



white dwarf PG 1351H-489 and found strong evidence of pe- 
riod doubling. Period doubling is a scenario where a sys- 
tem that primarily oscillates with a fundamental frequency 
vo, after change o f the contr o l parameter oscillates with fre - 



quency vo/2 (see iHelleman . , Il980l: iGurel & Rosslel Il979h . 



til the pioneering work of Baker et al. ( 19661) . which showed tems. ISerre et al 



Cascades of period doubUng is a typical transition routes to 
chaot ic behaviour . Such an effect is typical of chaotic sys- 
19961) analysed long time-series of photo- 



that a simple one-zone model of stellar oscillators with a 
nonlinear adiabatic term can in some specific conditions pro- 
du ce irregular pulsat i on. Similar results were also obta ined 
bvlMoore & Spiegel (Il966h . lAuvergne & Baglin I ( 119851) and 
Auvergne (Il988l) . The disadvantage of these models is their 



metric observations of the variable star R Scuti, one possible 
candidate pulsating star which exhibits a strange irregular light 
curve. By using the method of global flow rec onstruct i on (se e 
more about global time series predictions in iBrown I (1 1993b . 



Abarbanel et al. I (119941) ). they were able to reproduce the ob- 



dynamically unsta ble equil i brium state, whic h seems to be non 
realist ic for stars. Takeutil ( 1987 ) and later Saitou & Takeuti 



servation with the chaotic dynamics very well. Mo re details 



(Il989l) showed that similar, yet dynamically stable and more 



and ex amples of chaotic pulsations can be found in iBuchler 
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(Il987n . These clearly show that there is both theoretical and 
experimental evidence of quasiperiodic or chaotic behaviour in 
variable stars. The problem is how to distinguish between these 
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Fig. 1. Left panel: Chaotic attractor constructed from a numerical solution of the ODE system ([TJ with the parameters from Table 
[U - chaotic case. Right panel: Two-limit cycle attractor constructed from numerical solution of ([T]i with the parameters from 
Table[r| - periodic case. 



two possibilities, if only limited time series are available, often 
very unevenly sampled. In other words, what will be the result, 
if classical astronomical tools for the period searches, such as 
periodograms, the phase dispersion minimalization method or 
more sophisticated techniques are applied to data originating 
from quasiperiodic or chaotic sources. In the next paragraphs 
we tackle the problem using synthetic data. 

2. Generation of synthetic time series 

As an example of irr egular stellar variability d ue to pulsation 
we use the model of Tanaka & Takeutil ( 1988 ). In this model, 
the set of ordinary differential equations (ODE) for a stellar 
oscillator with nonlinear, nonadiabatic terms, reads 



df 
dy 
df 
dz 
df 



— — a + py + z 



= -jSy - pz- qy + syz 



(1) 



where the first two equations represent equations of motion and 
the last one is a state equation with a nonadiabatic term. We 
can adjust the behaviour of the model using different values of 
a,/3,p, p, q, s. The system of ODE of the first order was numer- 
ically solved with the Runge-Kutta method for two different 
sets of model parameters, the first one corresponding to the pe- 
riodic regime with a two-limit cycle attractor and the second 
one to the aperiodic, chaotic regime with a strange attractor 
similar to the Rossler type. The parameters are summarized in 
Table [U We use the initial conditions (xq - 0.37, y - -0.249, 
z = -0.174) for the chaotic solution. The initial conditions are 
very close to the orbit of the attractor in phase space. We ob- 
tain a time-dependence of the state variables x(f),y(f),z(f) in 
the form of a time series with an equidistant time step. We 
choose one state variable, namely x{t), as a representation of 
the measured variable. To simulate a synthetic signal from the 
quasiperiodic process, a different approach is employed. The 



signal is modeled by the following equation 

= fli sin (/if) + ai sin {fit) + a^ sin (/sf) (2) 

where a\,a2, 03 are amplitudes of incommensurate modes with 
frequencies /i,/2,/3. Their values are summarized in Tables [T] 
and|2] These parameters are carefully chosen in order to obtain 
similar dependencies of x(f) and Xquasi(?) for the chaotic and 
quasiperiodic case, namely amplitudes of oscillations (see Fig. 
|2]). The computed signal for each type of process is plotted in 
Fig. [2] Since we want to simulate real measurements, we add a 
synthetic signal with Gaussian noise. It is important that for a 
simulation of light curves from objects others than stars, such 
as QSOs, noise may by colored and then it is necessary to use 
a different type of noise, e.g. Poisson noise. The noise can be 
described using a normal distribution 



G{x) 



■V2^ 



exp((x - x)^/2cr^) 



(3) 



where cr is the variance and x is the mean value of the mea- 
surement errors. For numerical simulations of the noise with a 



normal distribution, we use the method of iHamminpl jl962l) . 
It was necessary to set the different values of cr for each type 
of signal in order to have similar S/N ratios. Data for synthetic 

Table 1. Parameters of the stellar pulsation model 
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P 
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cr 


X 


periodic 


-0.5 


0.5 


0.5 


3.2 


0.5 


1.0 


0.2 





chaotic 


-0.5 


0.5 


0.5 


4.0 


0.5 


1.0 


0.3 






time series satisfy the condition f > 300, i.e. reasonably far 
away from the starting point (f - 0). It is necessary to make 
sure that the solution is relaxed away from the influence of the 
initial conditions. 

It can be readily seen that chaotic and quasiperiodic oscil- 
lations (see Fig. O have very similar patterns, different to the 
simple periodic signal (two-limit cycle). 
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Fig. 2. Upper left panel: x{t) solution of the ODE system ([T]i for the periodic case with a two-limit cycle. Upper right panel: x(t) 
solution of equation Q for the quasiperiodic case. Lower left panel: x(t) solution of the ODE system ^ for the chaotic case. 
Lower right panel: Gaussian noise (|3]l with parameters o" = 0.1 and x -0. All physical quantities are in arbitrary units. 



3. Analysis of the signal 

3.1. Classical period analysis 

The most common techniques used for period searches are 
the power spectral analysis (PWS), the modified Lomb peri- 
odogram (LPD), the phase dispersion minimalization (PDM), 
and the string length method (SLM). In the first case, discrete 
Fourier transformation is applied to a d ataset resulting in a 



power spectrum (for more details see e.g. IScargle L 1 1 9821) . For 



unevenly sampled series it is better to us e a modificati on of 
PWS, method LPD (see for more details iLombL Il976h . 



The 



third method binds data in a grid of phases a nd minimizes 
phase dispersion statistics (se e for more details IStellingwerf , 
19781: lLafler& Kinman LI 1965b. The SLM method is similar to 



the previous one fsee iDworetsky L 1 1 983h . and the LPD method 



behaves like the PWS. So we limit on the two methods, PWS 



Table 2. Parameters for the quasiperiodic signal 



ai 



«3 



/i fl 



0.4 0.6 0.5 VI VI VI 0.1 



and PDM only. Obviously, PDM and PWS techniques applied 
on a periodic signal detect corresponding frequencies without 
any doubt, but, as we will show, a different situation arises for 
quasiperiodic and chaotic regimes. 

7. Quasiperiodic signal. The PWS analysis of the 
quasiperiodic signal produced by known incommensurate fre- 
quencies (see Table |2]i detects the frequencies used. Moreover, 
it is possible to recognize the leading frequency, which is in 
our case /2 = V3. The noise in the periodogram in Fig. |4] 
is caused by the simulated Gaussian noise of the signal. The 
PDM method yields a similar estimate of the input frequen- 
cies, and the results also show some other frequencies that are 
subharmonics or harmonics of the frequency set. The phase 
diagram in Fig.|5](on the left) is folded with the leading fre- 
quency /2. Clearly, the structure is not a simple sine wave. It is 
a quasiperiodic, and hence it is composed of a set of harmon- 
ics with slightly different phase. After accounting for Gaussian 
noise, this difference is smeared out. 

2. Chaotic signal. The results become much more interest- 
ing for the chaotic signal. Again, both PWS and PDM meth- 
ods detect periods, as can be clearly seen in Fig. H] After 
some checks we can recognize the two longest independent 
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Fig. 3. correlation integrals (|4|i for different lengthscales e and embedding dimensions m. The quasiperiodic case is in the left 
panel, chaotic case is in the right panel. 



periods. The power spectrum of the chaotic signal shows the 
most significant frequencies together with the frequency noise. 
Nevertheless, low-level noise is present even if Gaussian noise 
is not applied. It is the consequence of the chaotic behaviour 

Theoretically, power spectra of the chaotic data are com- 
posed of an infinite number of frequencies. This is true, how- 
ever, only if there is an infinite number of observations (see 
Regev L I2OO6I chap. 5). In reality, there is only a limited num- 
ber of observations, which are, furthermore, unevenly spaced. 
The phase diagram constructed for the most significant fre- 
quency derived /i = 0.1 179 is shown in Fig. |5](on the right). 
In the phase diagram without Gaussian noise, we can recog- 
nize a two-wave structure in the phase plot. On the other hand, 
when Gaussian noise is applied a simple noisy sine wave can 
be identified. This is similar to the quasiperiodic case results, 
although the data originated from a nonlinear chaotic process, 
where period has no meaning. Other frequencies fo,/^ lead to 
even less clear phase diagrams. 

The above results clearly show that it is very difficult to dis- 
tinguish between chaotic and quasiperiodic cases using period 
analysis only. Thus, we must look for other methods, which can 
help to resolve both processes. 



3.2. Nonlinear analysis 

For detailed analysis of the signal that can be generated by non- 
linear chaotic processes, we can use methods of nonlinear time 
series analysis. These methods were developed for classifica- 
tion and identification of chaotic ph enomenon in time seri es of 



detail iKantz & Schreiberl (|2004 ) as 



the experimental data (see review lAbarbanel et al. L 119931) . As 



a indication of chaotic behaviour we can use topological in- 
variants, namely the correlation dimension and Lyapunov ex- 
ponents. Algorithms for numerical e stimates of such invari- 
ants fr om time series are describe d by lAbarbanel et al. I (Il993h 



C(m, e) 



(A^ - «min)(A^ - n„ 



1) 



and bv lKantz & Schreiberl ( 120041) . In our case the indicator of 
chaotic behaviour is the correlation integral defined (see more 



(4) 

where m is the embedding dimension, Sn denotes the delay em- 
bedding vector, e is the lengthscale parameter , and is the 
number of points. A possible problem with temporal correla- 
tions is overcome by excluding pairs that are too close in time 
(parameter «min). Before we apply the procedure for calcula- 
tion of the correlation dimension to the noisy signal (chaotic 
or quasiperiodic), it is necessary to clean the noisy data. This 
is done by using nonlinear noise reduction with locally con- 
stant approximati ons - the simple nonlinear filter described by 
Schreiberl ( Il993h . The results of the calculation are shown in 
Fig. |3j The chaotic processes are characterized by a plateau 
in the plot of C(m, e) for different embedding dimensions (left 
panel Fig.O, which indicates a self-similar geometry. A linear 
fit to this plateau region leads to the estimation of correlation 
dimension 

Cdin, = 1.90 + 0.08. (5) 

This value agrees with the typical value of the Rossler type at- 
tractor In the case of the quasiperiodic signal we cannot iden- 
tify the plateau region (right panel Figl3ll. This indicates that 
there is no self-similar geometry for this type of the signal. 

4. Summary 

Analysis of synthetic signals simulating quasiperiodic and 
chaotic data led to the following findings: 

1) In the realistic case ( finite number of observations 
with Gaussian noise), it is impossible to distinguish between 
quasiperiodic and chaotic behaviour by classical methods of 
period analysis. 

2) When the data contain large amounts of noise, the 
chaotic or quasiperiodic character of the signal can be over- 
looked and a simple periodic solution with the strongest fre- 
quency can be determined from the period analysis. 
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3) In order to definitely rule out presence of chaos it is nec- 
cessary to use a nonlinear time series analysis. Thus, we must 
reconstruct the topologically similar phase portrait of t he sys - 
tem using the technique described by lAbarbanel et al. I ( Il993b . 
Then we must determine the invariants, such as the correla- 
tion dimension or Lyapunov coefficients from the reconstructed 
time series, and from their values we can decide if deterministic 
chaos is present in the data or not. However, both these proce- 
dures need better quality of data (noise, length of time series). 

4) We would like to draw attention to the fact that irregular 
chaotic behaviour may be more usual in astronomical data than 
thought before, and that some systems that were identified as 
quasiperiodic may actually be chaotic ones. 
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Fig. 4. Upper panels: Power spectra constructed for quasiperiodic oscillation (on the left) described by (|2]i and the chaotic case 
of stellar oscillation (on the right) described by ([TJ. Bottom panels: Periodogram constructed using the phase dispersion mini- 
malization method for the quasiperiodic case (pti the left) and chaotic case {on the right). For the frequency we used arbitrary 
units. 
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Fig. 5. Left upper panel: Phase diagram for the quasiperiodic case and the dominant frequency /i = V3 without noise. Right 
upper panel: Final phase diagram for the frequency f\ - 0.1176 of the chaotic case. Left bottom panel: The same as the left 
upper panel but with Gaussian noise. Right bottom panel: The same as the right upper panel, but with Gaussian noise. On the 
bottom panels both diagrams are similar, showing a simple sine-wave pattern. All physical quantities are in arbitrary units. 



